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ABSTRACT 


A method  has  been  developed  to  predict  the  lift  and  thrust  of  an  ejector 
wing  by  iterating  between  a viscous  solution  for  the  turbulent  entrain- 
ment of  the  primary  jets,  and  an  inviscid  solution  for  the  ejector  wing 
flow  field.  A two-dimensional  analysis,  which  utilizes  a turbulent 
kinetic  energy  model  for  the  jet  mixing  calculation  and  a higher  order 
panel  method  for  the  inviscid  flow  calculation  is  described.  The 
•complete  ejector  wing  geometry  is  analyzed.  Detailed  surface  pressures 
both  inside  and  outside  the  ejector  can  be  calculated.  A sample  cal- 
culation for  a typical  ejector  wing  configuration  is  compared  to  experi- 
mental .data. 


*The  work  reported  in  this  paper  was  supported  partly  by  NASA-Ames 
Contract  NAS2-10681,  ONR  Contract  N00014-78-C-0557,  and  partly  by 
the  Independent  Research  and  Development  Program  of  Rockwell  International 
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INTRODUCTION  • 


Although  analytic  methods  are  necessary  for  conceptual  studies  and  to 
reduce  test  requirements,  there  is  no  satisfactory  theory  for  predicting 
ejector  wing  performance.  Methods  have  been  developed  to  calculate  the 
surface  pressure  distributions  induced  by  an  ejector  of  given  thrust. 

These  calculations  are  based  on  the  no w classical  vortex  sheet  model  of 
the  pure  jet  flap  devised  by  Spence  (ref.  1).  Linearized,  thin  airfoil 
models  of  the  ejector  wing  were  developed  by  Chan  (ref.  2)  and  Woolard 
(ref.  3),  who  added  a sink  on  the  upper  surface  of  a jet-flapped  wing 
to  represent  the  entrainment  into  the  ejector.  Wilson  (ref.  4)  extended 
this  approach  by  including  the  effects  of  thickness  and  camber,  as  well 
as  deflection  of  the  jet  wake.  More  recently,  Dillenius  and  Mendenhall 
(ref.  5)  studied  three-dimensional  effects.  In  ail  these  methods-,  experi- 
mental data  are  used  to  specify  the  variation  in  ejector  thrust  during 
conversion  from  hover  to  conventional  flight. 

Such  an  empirical  approach  is  useful  for  estimating  pressure  distri- 
butions or  for  performing  parametric  analysis.  However,  a theory  for 
predicting  both  the  thrust  augmentation  ratio  and  the  initial  thrust 
angle  is  necessary  to  evaluate  significant  design  changes,  or  new  con- 
figurations for  which  there  is  no  data  base.  Bevilaqua  and  DeJoode 
(ref.  6)  developed  a method  for  predicting  the  thrust  augmentation  of 
stationary  ejectors  by  iterating  between  a viscous  solution  for  the 
entrainment  of  the  primary  jets  and  an  invi'scid  solution  for  the  pressures 
induced  on  the  ejector  duct  by  the  entrained  flow.  The  purpose  of  this 
paper  is  to  describe  cn  extension  of  this  analysis  developed  to  predict 
the  lift  and  thrust  of  the  more  complex  ejector  wing  in  transition  flight. 

Calculation  of  the  flow  field  is  more  difficult  in  this  case  because 
deflection  of  the  exhaust  jet  by  the  free  stream  influences  both  the 
thrust  of  the  ejector  and  the  lift  of  the  wing.  A two-dimensional 
analysis,  which  utilizes  a turbulence  kinetic  energy  model  for  the  jet 
mixing  calculation  and  a higher  order  panel  method  for  calculating  the 
surface  pressures,  has  been  developed. 

A description  of  this  analysis,  together  with  a computed  result 
is  given  in  the  following  sections. 


MATHEMATICAL  MODEL 


The  interactions  between  the  ejector  and  wing  flow  fields  are  corm 
puted  without  solving  the  full  Navier-Stokes  equations  by  iterating 
between  a viscous  solution  for  the  flow  through  the  ejector  and  an 
inviscid  solution  for  the  flow  around  the  wing.  Because  the  primary 
direction  of  flow  is  through  the  ejector,  the  governing  elliptic 
equations  can  be  reduced  to  a parabolic  set  which  can  be  solved  by 
marching  through  the  ejector  in  the  streamwise  direction.  The  mixing 
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of  the  turbulent  jets  can  then  be  used  to  define  an  equivalent  sink 
distribution.  The  requirement  that  the  ejector  wing  and  jet  boundaries 
must  be  streamlines  of  the  flow  together  with  appropriate  jet  dynamic 
boundary  conditions  determines  the  strength  of  the  wing  circulation. 

The  circulation  and  jet  shapes,  .in  turn,  control  the  entrainment  cal- 
culated in  the  next  iteration  for  the  viscous  solution. 


Potential -Flow  Calculation 

The  geometry  of  the  two-dimensional  ejector  wing  considered  in 
this  study  is  shown  in  Figure  1.  It  has  three  main  components:  a 
plain  flapped  wing,  a central  nozzle,  and  an  aft  flap.  Primary  jets 
are  injected  at  the  knee  of  the  forward  flap,  the  trailing  edge  of 
the  central  nozzle,  and  near  the  leading  edge  of  the  aft  flap.  The 
three  jets  grow  and  merge  to  define  the  jet  wake,  which  has  higher 
total  pressure  than  the  free  stream  does.  Consequently,  the  flow  con- 
sists of  two  regions  with  different  total  pressure,  thus  it  is  inhomo- 
geneous 

According  to  Kuchemann  and  Weber  (ref.  7),  the  inhomogeneous  flow 
may  be  made  homogeneous  without  changing  the  flow  velocity  field  by 
subtracting  the  total  pressure  difference  between  the  jet  wake  and 
.the  main  stream  from  the  static  pressure  inside  the  jet  wake.  In  the 
resulting  homogeneous  flow,  the  jet  boundaries  are  unknown  and  to  be 
determined  as  part  of  solution  subjected  to  the  usual  tangential  flow 
(kinematic)  condition  and  an  additional  dynamic  boundary  condition 

AH  = pUy  (1) 

where  p is  the  fluid  density  (both  jet  and  main  stream),  U is  the  mean 
velocity  of  the  jet  boundary,  y is  the  velocity  difference  across  the 
jet  boundary,  and  aH  is  the  total  pressure  difference  between  the  main 
stream  and  jet  wake. 

However,  if  the  jets  are  not  completely  mixed  by  the  ejector  exit, 
each  jet  before  complete  mixing  is  treated  as  a thin  jet  using  the 
classical  thin  jet  theory.  According  to  Spence  (ref.  1),  the  static 
pressure  jump,  Ap,  across  a thin  jet  is  balanced  by  the  rate  of  change 
of  jet  momentum,  J,  due  to  jet  curvature,  1/R,  and  is  related  to  the 
vortex  strength,  y,  (equivalent  to  the  velocity  difference  across  the 
jet)  of  the  jet  sheet  as  follows: 


Ap  = PUy  = -^  (2) 

The  solution  of  a thin  jet  is  obtained  by  satisfying  this  dynamic 
boundary  condition  together  with  its  kinematic  boundary  condition.  All 
the  flow  singularities  which  determine  the  potential -flow  solution  are 
shown  in  Figure  1.  The  cross  line  source  is  added  to  combine  with  the 
line  sink  to  simulate  the  doublet  effects  on  the  flow  field  for  modeling 
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the  effect  of  the  jet,  which  acts  like  an  elongated  actuator  disk  to 
draw  air  through  the  ejector. 

The  potential  flow  just  described  is  calculated  by  using  the  Hess 
(ref.  8)  higher-order  panel  method.  Both  airfoil  and  jet  boundaries 
are  defined  by  a series  of  discrete  points,  so-called  corner  points, 
as  shown  in  Figure  2.  Between  two  successive  corner  points,  the  true 
geometry  is  approximated  by  a curved  parabolic  panel.  A linear  vortex 
distribution  and  a linear  source  distribution  are  placed  on  each  of 
these  panels.  Source  singularity  strengths  are  chosen  a priori  and  set 
equal  to  twice  the  local  jet  entrainment  velocity.  Vortex  singularity 
strengths  and  the  jet  shapes  are  determined  by  satisfying  the  airfoil 
and  jet  kinematic  boundary  conditions  and  the  jet  dynamic  boundary  con- 
ditions, equations  (1)  and.  (2). 

Since  the  jet  dynamic  boundary  conditions,  equations  (1)  and  (2), 
are  nonlinear  and  the  jet  shapes  are  not  known,  a priori,  an  iterative 
procedure  shown  in  Figure  3 is  adopled  to  obtain  the  potential -flow 
solution.  Details  of  the  computational  procedure  are  given  in  ref.  9. 


Jet-Mixing  Calculation 

The  jet-mixing  calculation  is  a partially  parabolic  method  described 
in  detail  in  ref.  10.  The  flow  governing  equations  are  derived  from 
Reynolds'  equations  for  turbulent  flows,  by  neglecting  streamwise  dif- 
fusion and  including  curvature  effects.  A TKE  turbulent  model  (ref.  11) 
modified  to  include  the  curvature  effect  is  used  to  determine  the 
turbulent  viscosity.  A set  of  finite-difference  equations  are  formed 
by  integrating  the  governing  differential  equations  over  a small  control 
volume.  The  resulting  finite  difference  equations  are  solved  iteratively 
for  velocity  and  pressure  fields.  Briefly.,  the  iterative  procedure  begins 
with  an  initial  guess  of  the  pressure  field,  solves  the  momentum  equations 
for  the  velocity  components  using  a triadi agonal  matrix  algorithm,  corrects 
the  pressure  and  velocity  fields  to  satisfy  the  continuity  equation,  and 
repeats  the  process  until  convergence  is  obtained. 


SOLUTION  MATCHING  PROCEDURE 


The  solution  is  iterative,  since  the  two  individual  flow  problems 
and  their  coupling  are  nonlinear.  The  iterative  procedure  used  is 
summarized  schematically  by  the  flow  chart  of  Figure  4.  Details  of  the 
computational  procedure  are  described  in  ref.  9.  Presently,  the  potential - 
flow  program  and  the  jet-mixing  program  are  separated  and  communicate 
through  external  data  transformation. 


EXAMPLE  SOLUTION 


To  test  the  present  matching  solution  procedure  for  ejedtor  wings, 
the  model  wing  "Configuration  Augmenter  1"  (Figure  7)  of  ref.  12  was 
analyzed  at  a transition  operating  condition.  The  tunnel  velocity  was 
34.4  m/sec  (113  ft/sec)  for  a tunnel  static  pressure  of  1.024  atmospheres. 
The  effective  angle  of  attack  was  2 degrees  and  the  momentum  coefficient 
was  2.  Figures  8 through  12  show  the  calculated  results  and  Stewart's 
(ref.  12)  experimental  data.* 

Comparison  of  the  distributions  reveals  that  the  differences  between 
the  experimental  and  theoretical  pressure  are  dramatically  reduced  when 
•jet  entrainment  effects  are  included  in  the  calculations.  Examination  of 
the  experimental  data  shows  large  discrepancies  near  tie  leading  edge  and 
on  the  upper  surface  of  the  centerbody.  These  are  most  likely  due  to 
flow  separation  in  these  regions.  The  flow  separation  may  be  due  to 
the  high  local  angle  of  attack  as  indicated  by  the  calculated  flow 
streamline  pattern  about  the  ejector  wing  shown  in  Figure  12.  The 
present  computer  program  cannot  calculate  a separated  flow. 


CONCLUSIONS 


1.  A viscid/inviscid  interaction  analysis  has  been  developed  to 
predict  the  thrust  augmentation  ratio  and  initial  thrust  angle  of 
ejector  wing  configurations.  This  provides  an  advance  over  classical 
methods  of  analysis,  which  require  these  parameters  to  be  specified 
as  an  input. 

2.  Comparison  of  the  predicted  surface  pressure  distributions  with 
experimental  data  establishes  confidence  in  the  model  of  the  flow 
field.  But  in  addition,  a greater  understanding  of  the  performance 
of  ejector  wings  has  been  obtained  from  the  analysis. 

3.  Examination  of  the  computed  streamlines  ana  surface  pressures  reveals 
the  somewhat  surprising  result  that  the  forward  stagnation  point  is 
located  near  the  trailing  edge  of  the  forward  flap.  This  is  a 
result  of  the  large  circulation  induced  by  the  jet  flap  effect 

and  suggests  that  a leading  edge  device  may  be  required  to  achieve 
maximum  lift  and  thrust.  Further,  the  surface  pressure  distributions 
on  the  other  two  elements  suggest  that  the  flow  is  more  likely  to 
separate  from  the  central  nozzle  than  the  aft  flap.  Since  the 
flow  into  the  ejector  is  being  accelerated  past  the  nozzle,  this 
was  also  an  unexpected  result. 

4.  The  present  analysis  can  be  improved  by  including  boundary-layer 
and  three-dimensional  effects. 

5.  It  would  be  useful  to  have  test  data  for  a two-dimensional  configuration,  ^ 
since  no  comparison  with  available  three-dimensional  data  can  be  exact!  < 
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